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ABSTRACT 

We present the first observations of extragalactic pulsating stars in the K2 ecliptic survey of the 
Kepler space telescope. Variability of all three RR Lyrae stars in the dwarf spheroidal galaxy Leo IV 
were successfully detected, at a brightness of Kp ^21.5 mag, from data collected during Campaign 1. 

We identified one modulated star and another likely Blazhko candidate with periods of 29.8±0.9 d and 
more than 80 d, respectively. EPIC 210282473 represents the first star beyond the Magellanic Clouds 
for which the Blazhko period and cycle-to-cycle variations in the modulation were unambiguously 
measured.The photometric [Fe/H] indices of the stars agree with earlier results that Leo IV is a very 
metal-poor galaxy. Two out of three stars blend with brighter background galaxies in the K2 frames. 

We demonstrate that image subtraction can be reliably used to extract photometry from faint confused 
sources that will be crucial not only for the K2 mission but for future space photometric missions as 
well. 

Subject headings: stars: variables: RR Lyrae — methods: observational — techniques: photometric 


1. INTRODUCTION 

Kepler has provided incredible results on extrasolar 
planets and planetary systems, as well as on stellar as¬ 
trophysics. The space telescope was designed to be the 
most precise photometer ever built in order to detect 


the t ransits of numerous small planets (Borucki et al. 
2010). To achieve this goal, the original mission of Ke¬ 


pler focused almost exclusively on the Galactic stellar 
population, approximately 170 000 targets, within the 
Lyra-Cygnus field. But after the failure of two reaction 
wheels the telescope permanently lost its ability to main¬ 
tain its original attitude. Soon a new mission, called K2, 
was initiated to save the_othe rwise healthy and capable 
space telescope (Howell et al.|2014). Since then, Kepler 
has been observing m shorter, 75 day long campaigns 
along the Ecliptic to balance the radiation pressure from 
the Sun. 

The new fields opened the possibility to extend the 
capabilities of Kepler into new are as an d we set out to 
explore the limits of the K2 mission. Molnar et al. (2015) 
investigated the fir st obser vations of field R R Lyrae stars. 


Szabo et al. (2015) and Pal et al. (2015) showed that 
Kepler can be used for Solar System research, including 
the detection of main-belt asteroids and trans-neptunian 
objects (TNOs). In this paper we continue this work, ex- 
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panding the reach of Kepler towards extragalactic pul¬ 
sating stars. The three RR Lyrae stars in the galaxy 
Leo IV are the first non-cataclysmic stellar targets Ke¬ 
pler ever detected outside the Galaxy. We note in passing 
that the space telescope has already o bserved fo ur super¬ 
novae during the original mission (Oiling et al.||2015). 

Leo IV is one of the small and faint dwarf spheroidal 
galaxies around the Milky Way that were recently dis¬ 
covered with the help of the Sloan Digital Sky Survey, 
with a mass and luminosit y of M = ( 1.4 ± 1. 5) • IQ 6 M^ 
and My = —5.1 ± 0.6 mag 


(Belokurov et al.|2007l jSimon 
ted at a heliocentric distance 


& Geha 2007). It is located" 
of 154 ± 5 kpc, and has a half-light radiu s of 3.3 arcmin 


translating to a physical size of 160 pc (Moretti et al. 


2009). Most of the stars in Leo IV are older than 12 Gyr. 


although another, brie f epis ode of star formation likely 
occurred 1-2 Gyr ago (Sand et al.|[2010). The galaxy is 
very metal-poor: different studies determined th e aver 


& Geha 2007 

Sand et al. 2010 

) or possibly as low as 

\tr/H\r= -2 

( 

Kirby et al. 

2008). The galaxy was 

also surveyed by M 

.oretti et al. 

(2009) who discovered 


four variable stars, three fundamental-mode RR Lyrae 
(RRab) and a single SX Phe pulsator. 

RR Lyrae stars are ubiquitous in all nearby galaxies 
and have been detec ted even beyond the Local Group 
(Da Costa et al. 2010), but our knowledge decreases with 
distance. Our closest neighbors, the Magellanic Clouds 
and the Sagittarius dwarf and its stellar stream have been 
surveyed extensively both in terms of stars and tempo- 


grams (see, for example, 

Alcock et al. 

1996 

Soszynski 

et al. 

2009 2010 

2014). 

But studies oj 

t other galaxies, 


amounts of data. These observations are usually aimed 
at distance determination and population studies. 

Extragalactic surveys have their limitations, how¬ 
ever. Mu ltiple studies revealed several candidate Blazhko 


stars. 

24 stars 


Stetson et al. (2014|) for example, identified 
with signs of amplitude variation out of 194 
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Fig. 1.— The K2 target pixel masks covering parts of Leo IV. 
The three RR Lyrae stars are indicated. The insert shows high- 
resolution deep imagery obtained by combining data from the 
William Herschel Telescope, the Isaac Newton Telescope, andjthe 
Southern Astrophysical Research telescope, collected by |Moretti et| 
|al.|(|2009). Red circles repre^sent the apertures used in the reduction 
ot the JA 2 data (see Figure [ 2 ). 


RR Lyrae s tars in Leo I, and we know one candidate 
in M31 too (Br own et al.| |2007). These low occurrences 
are likely just lower limits: we now know that ab out 50% 
of RRab stars are modulated in the Milky Way (Jurcsik 
et al. 2009 Benko et al. 2014). But in order to deter- 
mme the accurate occurrence rates and modulation peri¬ 
ods of Blazhko stars, one usually requires several weeks 
of intensive observations, and thus unambiguous detec¬ 
tions have not yet been made. Also, data from Kepler 
revealed that continuous observations are the key to de¬ 
tect any cycle-to-cycle variations—ano ther feature that 
earlier studies lacked (Szabo et al. 2010). Therefore we 
set out to gather the first extended, uninterrupted obser¬ 
vations from extragalactic pulsating stars with the help 
of K2. 


2 . OBSERVATIONS AND DATA REDUCTION 

Leo IV fell into the field-of-view of the first scientific 
campaign (Cl) of the K 2 mission and was observed be¬ 
tween May 30 and August 21 , 2014, corresponding to 
Barycentric Julian Day (BJD) 2456808.18-2456890.33. 
We requested to observe the 3 RR Lyrae variables and 
the 16 brightest giant stars in the galaxy through the K2 
Guest Observer proposal GO 1019 in long cadence mode, 
with an integration time of 29.4 min. We estimated the 
maximum brightness of the RR Lyrae stars to be around 
Kp = 21 — 21.2 mag in the Kepler passband, leading to 
a precision of a few tenths of a magnitude per long ca¬ 
dence data point. At that time the much fainter SX Phe 
star (V = 23.0 mag) seemed to be beyond the capabil¬ 
ities of Kepler so we did not include it in the proposal, 
although the recent experiences with the similarly faint 
TNO, 2002 G V 31 , suggest that we could have detected its 
variation too ( Pal et al.|2015 ). In this paper we focus on 


TABLE 1 

Summary of the K2 observations of the RR Ly targets in 

THE DWARF GALAXY LEO IV. 


Object 

Mask 

(px) 

R.A. 

(J2000) 

Decl. 

(J2000) 

V 

(mag) 

ID 

210282472 

15x14 

ll h 32 m 55.8 3 

—0°33'29.4" 

21.46 

V2 

210282473 

15x13 

ll h 32 m 59.2 s 

—0°34 , 03.6" 

21.47 

VI 

210282474 

14x14 

ll h 33 m 36.6 s 

—0°38'43.3" 

21.52 

V3 


Note. — IDs refer to the identifications given by |Moretti et al.| 
( |2QQ9| ). 


the three RR Lyrae stars from the sample. The summary 
of their observation parameters is shown in Table [l] 

For all of the tasks described below we employe d excl u- 
sively the utilities shipped within the Fitsh 8 (|Pal|2012). 
Fitsh is a lightweight, yet comprehensive, fully open 
source astronomical data reduction and analysis software 
package. It is comprised of a collection of standalone bi¬ 
nary programs that are utilized through various UNIX 
shell scripts. The accurate photometric time series for 
these three extragalactic Leo IV RR Lyr targets have 
been obtained as follows. The image scale of Kepler is 
3.98"/px, and the point spread functions are at least 3- 
5 px wide that leads to confusion of nearby sources. Al¬ 
though Leo IV is sparse enough to be resolved even with 
Kepler , two out of the three RR Lyr targets were ex¬ 
tremely close to bright background galaxies (i.e., nearly 
within a pixel, see Figures [l] and |2|, accurate photometry 
can only be done by involving image subtraction tech¬ 
niques. Since the pointing jitter of Kepler in the K2 mis¬ 
sion was in the range of a pixel, the individual frames had 
to be adjusted to the same reference system before per¬ 
forming any kind of differential image analysis. In order 
to precisely estimate the shifts and rotations between the 
(subsequent) frames, we involved additional K2 frames 
since the raw K2 data do not include a frame-by-frame 
centroid coordinate. These additional 14 fields were the 
target pixel files of EPIC 201424914, 201430029 and the 
12 stamps having an identifier between 210282475 and 
210282486, all from the GO1019 proposal (where EPIC 
refers to the K2 Ecliptic Plane Input Catalog). These 
3 + 14 stamps were combined to a single image and the 
foreground stars were used to obtain precisely the trans¬ 
formation between these images (Figure [l]). 

Since the K2 mission observes along the ecliptic plane 
in an approximately 10°-wide area, the background of 
these observations also varies gradually due to the in¬ 
creasing amount of zodiacal light throughout the cam¬ 
paign. Hence, prior any differential analysis, the back¬ 
ground must also be subtracted. Luckily, the fields con¬ 
taining the 3 RR Lyr targets as well as the additional 
14 stamps were not crowded, and the background can 
simply be determined by considering the median of all 
observed pixel values. 

After the frames were registered to the same reference 
system and the background was subtracted accordingly, 
we chose every 20th image to obtain median-averaged 
master frame used as a reference for image subtraction. 
Since the number of stars having a good signal-to-noise 
ratio in these frames were insufficient and the instrument 

8 http://fitsh.szofi.net/ 






































3 


Extragalactic RR Lyrae stars with K2 



Fig. 2.— Image stamps showing the vicinity of the K2 extragalac¬ 
tic RR Lyrae targets 210282472, 210282473 and 210282474. The 
first row shows the stamps on the master median-combined image 
used in the process of differential photometry. The second row 
shows stamps created as the average of a roughly dozen differen¬ 
tial images taken at the brightest phase of the RR Lyr oscillations. 
The photometric aperture with an r = 1.5 px radius used in the 
procedure is shown as a red circle. Stamps in the third row show 
the respective SDSS DR9 images. All of these stamps cover an 
area of 64" x 64" on the sky, equivalent to 16 X 16 Kepler pixels. 


does not show any significant variation in its point-spread 
function, we did not perform any cross-convolution be¬ 
tween the adjusted images. This procedure also simpli¬ 
fied the evaluation of the photometry and hence it was 
not ne cessary to perform photometry on convolved aper¬ 
tures (Pal 2009). Therefore, the final fluxes were ob¬ 
tained as the sum of the flux in the subtracted frame 
and in the master frame, i.e., Fi = + Fmaster- 

For these flux estimations we used simple aperture pho¬ 
tometry with an aperture radius of r = 1.5 px (Figure 
[2|. However, we must note that the reference flux on the 
master frame can only be estimated with an additional 
systematic uncertainty in the case of the confused sources 
(210282472 and 210282473). Namely, we compared the 
peak pixel values on the stellar sources with the fluxes 
obtained using an aperture of r = 1.5 px. The ratio of 
these two numbers aid the estimation of the reference 
flux in the case of the confused sources as well and we 
found that the systematic uncertainty of this method is 
in the range of 2-3% for these sources. Due to the nature 
of the differential photometry, this kind of systematics 
yields a systematic amplification in the variability ampli¬ 
tudes. The random errors in the photometric fluxes have 
been estimated using the photon noise derived from the 
instrumental gain of 113.3e _ /ADU as well as the stan¬ 
dard deviation of the background pixels near the target 
sources. Due to the differential processing, this latter 
component also includes the photon noise of the nearby 
confusing sources. 


TABLE 2 

K2 PHOTOMETRIC DATA OF THE RR LYR STARS IN LEO IV. 


Object 

Time 

(BJD) 

Brightness 
(.R mag) a 

Error 

(mag) 

Type 

210282472 

2456810.28292 

20.865 

0.062 

raw 

210282472 

2456810.30336 

20.659 

0.048 

raw 

210282472 

2456810.32379 

20.974 

0.058 

raw 

210282473 

2456810.28292 

21.247 

0.069 

raw 

210282473 

2456810.30336 

21.188 

0.069 

raw 

210282473 

2456810.32379 

21.240 

0.076 

raw 

210282474 

2456810.28292 

21.821 

0.122 

raw 

210282474 

2456810.30336 

21.585 

0.093 

raw 

210282474 

2456810.32379 

21.565 

0.085 

raw 

210282472 

2456810.28292 

20.886 

0.062 

proc 

210282472 

2456810.32379 

20.995 

0.058 

proc 

210282472 

2456810.34422 

21.004 

0.054 

proc 


Note. — Table [2] is published in its entirety in the elec¬ 
tronic edition of the Astrophysical Journal Letters. A portion 
is shown here for guidance regarding its form and content. 


a Magnitudes shown here are transformed to USNO-Bl.O R 
system, see text for further details. Indices “raw” and “proc” 
correspond to raw and processed (sigma-clipped and Fourier- 
filtered) data, respectively. 


3. LIGHT CURVE ANALYSIS 

Once the raw light curves were obtained, we removed 
the outliers from the light curves with a sigma-clipping 
algorithm. We subtracted an initial Fourier fit from the 
light curves, clipped the residual at the 3cr levels three 
consecutive times and then restored the original signal. 

Finally, we applied Fourier filtering to remove the slow 
variations. The photometry of these faint targets is very 
sensitive to the variations in the background flux level. 
Even the best raw light curves we extracted contained 
additional variations up to 0.1 — 0.2 mag in the last 
10 — 20 d of data. However, a simple high-pass frequency 
filter would remove any modulation signal along with 
the low-frequency noise. So instead we removed all low- 
frequency components below 0.5 d -1 that were stronger 
than 10 mmag and were not connected to any potential 
modulation frequencies. The final light curves are dis¬ 
played in the left panels of Figure [3] The right panels 
show the folded light curves, along with the binned phase 
curves. 

The increasing background from the zodiacal light low¬ 
ered the photometric accuracy of individual data points 
towards the end of the campaign. The accuracy of a 
single LC point was similar for EPIC 210282472 and 
210282473: 0.050-0.055 mag at the beginning and 0.078- 
0.092 mag during the last days of the measurements. 
The errors are somewhat higher for the unblended star 
EPIC 210282474, starting from 0.087 mag and reaching 
0.107 mag at the end of the campaign. 210282472 and 
210282473 also appear to be brighter than 210282474 
by 0.50 and 0.34 mag, respectively. As mentioned, the 
main component of this difference comes from the un¬ 
certainty of the zero-point determination for the con¬ 
fused sources. Similar differences appear in the peak- 
to-peak amplitudes that are lowered to 0.4 and 0.6 mag 
for 210282472 and 2102 82473, compa red to 0.94 mag for 
210282474. In contrast, Moretti et al. (2009) found the 
amplitudes to be quite similar: 0.64, 0 .y 3^ and 0.65 mag 
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Fig. 3.— Light curves of the three RR Lyrae stars. Left panels: light curves after outlier r emoval and Fourier filtering. Note the variable 
amplitude of EPIC 210282473 in the middle panel. Designations from Moretti et al. 
folded light curves (pink dots and lines) and the binned phase curves (black points). 


d. (|2009|) are also indicated. The right panels show the 
JT We used 75 bins per pulsation period for each star. 


in V band, respectively. All these differences between 
EPIC 210282474 and the other two stars indicate that 
the confusing sources contribute considerably to the mea¬ 
sured flux levels of EPIC 210282472 and 210282473 in¬ 
creasing their brightnesses and lowering the amplitudes 
and photon noise levels we measure. Nevertheless, this 
zero-point uncertainty only changes the overall ampli¬ 
tudes but does not distort the shape of the light curves, 
therefore has no further effect on the astrophysical con¬ 
tent of the data. 

For all of these three objects, photometric magni- 
tn des have be en transformed into USNO-Bl.O R system 
(|Monet et al.||2003|). This pr ocedure has been performed 
similarly as it was done by |Pal et al.| (|2015|) . The raw 
and processed photometric data series are displayed in 
Table [2] Note that the photometric errors shown in this 
Table do not include the aforementioned zero-point un¬ 
certainty. 

We carried out a s tandard Fourier ana lysis with the 
Period04 software (Lenz & Breger 2005), using multi¬ 
frequency least-squares fits and consecutive prewhiten¬ 
ings. The identified frequency components are listed 
in Table [3j We detected modulation triplets (n/i har¬ 
monics and corresponding nf\ db f m sidepeaks, where 
n = 1,2,3...) in the Fourier spectra of EPIC 210282473 
and 210282474. We note, however, that in the latter 
case the length of the data is insufficient to properly re¬ 
solve the side-peaks, and their distances correspond to 
the length of the data instead of the modulation period, 
so we did not include them in the fit. 

We could also identify the characteristic signal of the 
spacecraft attitude correction maneuvers, at f corr = 
4.08 d -1 and its harmonics, but with signal-to-noise ra¬ 
tios below 4. Apparently, very faint objects suffer less 
from the attitude changes of K2 than bright targets. 





Fig. 4.— Fourier transform of the light curve of EPIC 210282473. 
The top panel shows the original spectrum. Significant modulation 
sidepeaks are visible in the middle panel after prewhitening with 
the main peak and its harmonics. The insert shows the triplet 
around the position of the main peak (dashed line). The bottom 
panel shows the residual after the side peaks have been prewhitened 
as well. The red solid and blue dashed lines are the 3 and 4 
SNR levels, respectively. The fcorr label marks the position of 
the marginally detected spacecraft attitude correction frequency. 
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0.2 0.4 0.6 0.8 1 1.2 1.4 

Pulsation phase (P=0.61903 d) 


0.2 0.4 0.6 0.8 1 1.2 1.4 

Modulation phase (P=29.8 d) 


Fig. 5.— Blazhko effect in EPIC 210282473. Left: data folded 
with the pulsation period. Binned data from 2-day long sections 
are overlaid, red dots show the maximum-amplitude modulation 
phase, blue circles show the minimum-amplitude phase. Right: 
light curve folded with the modulation period. 


Other significant frequency peaks, such as low-amplitude 
additional modes or half-integer peaks related to pe¬ 
riod doubling—both observed in many modulated RRab 
stars—did not appear after prewhitening. This is not 
surprising, however, considering that the 3 SNR (signal- 
to-noise ratio) significance limit in the Fourier spectra is 
around 8-12 mmag, well above the strongest additional 
peaks detected in the original K epler sa m ple. The ad¬ 
ditional modes described by|Benko et al. (2014) reached 
only 3 mmag, so similar modes would be undetectable in 
the Leo IV stars. The Fourier spectra of EPIC 210282473 
are shown in Figure [4j 

Such faint targets are very sensitive to any instrumen¬ 
tal problems arising both from the observations or the 
data reduction. Counterintuitively, the light curve of the 
only unblended star, EPIC 210282474 was the most prob¬ 
lematic in our case. Even after the Fourier filtering, some 
slow, low-amplitude 0.1 mag) variations still persisted 
in the last 15 days of the light curve (Figure [3]). Exam¬ 
ination of the K2 full-frame images did not reveal any 
obvious contaminating sources in the same or adjacent 
modules. The SDSS images did not reveal any close-by 
objects either (see Figure [2|. We visually inspected the 
stacked image obtained using the observations taken with 
the 2.5 m Isaac Newton Telescope, the 4.2 m William 
Herschel Telescope, and the 4.1 m Southern Astro physi- 
cal Re search Telescope (SOAR) in 2007 by Moretti et al. 
7 2009), but found no significant nearby objects to within 
0". This image is displayed in Figure IT] Images of the 


other two targets indicate that the confusion is caused 
almost exclusively by a single bright galaxy near each 
star. Another star can be detected 2.1" from 210282472, 
but it is about 3.5 mag fainter than the RR Lyrae target 
so its contribution is negligible. 

4. RESULTS 

4.1. Blazhko effect in Leo IV 

One star, EPIC 210282473, is clearly modulated. Al¬ 
though the signs of amplitude changes in RR Lyrae stars 
have b een observed as far as the Andromeda galaxy 
(Brown et al. 2007), this is the first star beyond the 
Magellanic Clouds for which the detailed modulation 
properties can be determined. The K2 data covers al¬ 
most three Blazhko cycles. The distance of the triplet 
side peaks (n/i db f m ) suggested a modulation period of 
Pm = 30 zb 3 d, but the independent fits to the mod¬ 
ulation peaks resulted in large uncertainties. We also 
calculated the amplitude and phase variations (Ai, (f \) 


of the main frequency peak (/i) by dividing the data into 
short segments. A sine fit to Ai(£) itself led to a more 
precise Blazhko period: P m = 29.8 db 0.9 d. 

Figure [5] shows the modulation properties of 
EPIC 210282473. The right panel is the light curve, 
folded with the modulation period: the presence of the 
Blazhko effect is even more clear here than in the nor¬ 
mal light curve. The left panel illustrates how the shape 
of the light curve itself changes over a modulation cy¬ 
cle. We selected two 2-day long sections of a low- and 
a high-amplitude state, respectively, and calculated the 
phase-binned values of the data. Red squares and blue 
circles show the two states, respectively: the differences 
in amplitude are evident around maximum and minimum 
light, but the phases at which the extrema occur show 
little to no shifts. 

Another representation is included in Figure [6] where 
the upper and middle panels show the variation oF the A\ 
and i Fourier-ter ms. These figures confirm that there is 
significant amplitude variation in the data, but the phase 
modulation is almost negligible in the star. Moreover, 
the middle column of Figure [6] indicates that the cycles 
are not repetitive and the amplitude and phase variation 
are not strictly correlated to each other, similar to some 
modulated stars in the original Kepler sample (Benko et 
aO2014|). 

The case of EPIC 210282474 is more ambiguous, since 
the light curve of the star suffers from instrumental ef¬ 
fects more than those of the other two. The Fourier- 
filtered light curve shows some amplitude and phase 
variation but those are suspiciously—but not exactly— 
symmetric to the mid-campaign data download period. 
This symmetry in the amplitudes invites to suspect some 
kind of instrumental origin, although we could not iden¬ 
tify any contaminating source. We also detect weak 
phase variations that seem to follow the changes in pul¬ 
sation amplitude (right column of Figure [6]). If we accept 
that these variations indeed originate from modulation, 
its period is clearly longer than the length of the cam¬ 
paign. We identified side-peaks in the frequency spec¬ 
trum, but given the limited length of the data, they are 
not resolved properly. Overall, since we detect the vari¬ 
ations in the first, better-quality part of the light curve 
as well, we accept this signal as a likely (although not 
unquestionable) detection of the Blazhko effect. Most 
stars in the Milky Way h ave m o dulat ion periods longer 
than 10-20 d (Szczygiel & Fabrycky 2007): the values 
we determined for the Leo IV stars agree with that pe¬ 
riod distribution. Although the Leo IV sample is very 
limited, it confirms that the Blazhko effect is abundant 
in other galaxies as well. 

EPIC 210282472 also shows some phase variations but 
only marginal changes in amplitude (left column of Fig¬ 
ure [6|. However, these amplitude changes can be partly 
attributed to the large scatter of the light curve points 
that make the phase determination uncertain for shorter 
data segments (10 days per segment in this case). Given 
the lack of significant amplitude variations, we classify 
this star as non-modulated. 

We also calculated the Fou rier p hase v alues for the 
original observations of Moretti et al. (2009) with the fre¬ 
quencies determined from the K2 data, and transformed 
them to the more classical O-C diagrams. The results, 
shown in the bottom row of Figure [6j suggest that pul- 
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TABLE 3 

Frequency tables for EPIC 210282472, 210282473 and 210282474. 


EPIC 

Freq. ID 

Freq (d 1 ) 

Amp (mag) 

0 (rad/27r) 

±F (d- 1 ) 

±A (mag) 

±cf) (rad/27r) 

210282472 

fo 

1.40872 

0.1400 

0.1786 

0.00011 

0.0023 

0.0024 

- 

2/o 

2.81745 

0.0640 

0.7512 

0.00023 

0.0023 

0.0053 

- 

3/o 

4.22617 

0.0427 

0.380 

0.00035 

0.0023 

0.008 

- 

4/o 

5.6349 

0.0217 

0.999 

0.0007 

0.0020 

0.016 

- 

5/o 

7.0436 

0.0106 

0.679 

0.0014 

0.0019 

0.032 

210282473 

fo 

1.61545 

0.2249 

0.1459 

0.00008 

0.0032 

0.0017 

- 

2/o 

3.23090 

0.0869 

0.6798 

0.00020 

0.0032 

0.0045 

- 

3/o 

4.84635 

0.0616 

0.2126 

0.00028 

0.0031 

0.0064 

- 

4/o 

6.46181 

0.0405 

0.801 

0.00042 

0.0031 

0.010 

- 

5/o 

8.0773 

0.0265 

0.349 

0.0006 

0.0030 

0.015 

- 

6/o 

9.6927 

0.0152 

0.950 

0.0011 

0.0028 

0.026 

- 

fo ~ fm 

1.5818 

0.0261 

0.016 

0.0006 

0.0025 

0.015 

- 

fo + fm 

1.6491 

0.0150 

0.344 

0.0012 

0.0031 

0.029 

- 

2/o fm 

3.1973 

0.0199 

0.485 

0.0008 

0.0030 

0.019 

- 

3/o fm 

4.8127 

0.0132 

0.998 

0.0012 

0.0030 

0.028 

- 

4/o - fm 

6.4282 

0.0104 

0.507 

0.0016 

0.0031 

0.037 

- 

5/o — fm 

8.0436 

0.0119 

0.144 

0.0014 

0.0031 

0.033 

- 

5/b + fm 

8.1109 

0.0098 

0.436 

0.0018 

0.0029 

0.041 

- 

6/0 + fm 

9.7263 

0.0095 

0.097 

0.0018 

0.0029 

0.043 

210282474 

fo 

1.58582 

0.2899 

0.7674 

0.00007 

0.0032 

0.0017 

- 

2/o 

3.17164 

0.1219 

0.9159 

0.00018 

0.0032 

0.0041 

- 

3/o 

4.75747 

0.0846 

0.0845 

0.00025 

0.0036 

0.0058 

- 

4/o 

6.34329 

0.0553 

0.300 

0.00039 

0.0035 

0.009 

- 

5/o 

7.9291 

0.0313 

0.516 

0.0007 

0.0035 

0.016 

- 

6/o 

9.5149 

0.0136 

0.677 

0.0016 

0.0035 

0.036 


Note. — Frequencies, amplitudes, phases and their respective uncertainties. The modulation frequency for EPIC 210282473 was set to 
fm = 0.03361 d _1 . 
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Fig. 6.— Top row: variation Ai, the amplitude of the fi frequency component in all three stars. Middle row: O-C variati on, calc ulated 
as the variation of the Fourier-parameter. Bottom row: long-term period drifts of all three stars between 2007 (Moretti et al. 2009 ) and 
2014, the K2 observations. 
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TABLE 4 

Pulsation periods, modulation periods and photometric 
[Fe/H] indices of the RR Lyr targets. 


Object 

Pulsation per. (d) 

Mod. per. (d) 

[Fe/H] 

472 

0.709862 ± 5.5 ■ 10“ 5 

— 

-2.16 ±0.10 

473 

0.619022 ±3.1 ■ 10“ 5 

29.8 ±0.9 

-2.64 ±0.10 

474 

0.630588 ± 4.4 • 10“ 5 

> 80 

-2.62 ±0.15 


sation period of EPIC 210282473 got shorter during the 
7.1 yr separating the two data sets while the periods of 
EPIC 210282472 and 210282474 got longer. Other ob¬ 
servations between the two epoch were too sparse to de¬ 
termine the pulsation phase and O-C values. Therefore 
actual rates of period change cannot be determined yet. 

4.2. Photometric metallicities 
A great advantage of RR Lyrae stars is that we can es; 


timate their [Fe/H] in di ces from photometry alone (Ju- 


rcsik & Kovacs 1996). Nemec et al. (2013) have call 


brated the relation between the light curve parameters 
and the spectroscopically determined metallicity values 
for the Kepler passband and showed that the photomet¬ 
ric and spectroscopic [Fe/H] indices agree to ±0.1 dex 
for all but the extremely modulated stars. We calcu¬ 
lated the photometric [Fe/H] values for all three stars 
with this method: the results are summarized in Table 
[4] The accuracy intrinsic to Fourier parameter determi¬ 


nation is ±0.03 dex. However, according to |Nemec et al, 
(2013), the uncertainties of the relation are dominated 


by" the calibration of the equation and the determina¬ 
tion of the spectroscopic values, so we assumed higher 
uncertainties for all stars (Table [4j 
Of the three stars, only EPKJ 210282472 has a spec 


troscopicall y determined index, [Fe/H] = —2.03 (Kirby 
et al. 2008) that is higher, but within the respective er¬ 


ror bars compared to our photometric ind ex of —2.16 ± 
0T0, assuming the same uncertainties. [M oretti et al.| 
(2009) determined the photometric metallicity for EPIC 
210282473, but their result, [Fe/H] = —2.11 is signifi¬ 
cantly higher than the K2-derived value (—2.64 ± 0.10). 
This discrepancy can lar gely be attributed to th e fact 
th at the linear relation of Jurcsik & Kovacs| (|1996|) , used 
by Moretti et al. (2009), overestimates the [Fe/Hj index 
by as much as ~ 0.3 dex for very metal-poor stars. The 
K2 data also revealed that the star is modulated, and 
the poor modulation phase coverage may have affected 
the original metallicity determination. Since the Blazhko 
effect in the star is well covered by the K2 data and does 
not show extreme amplitude changes or cycle-to-cycle 
variability, we cons ider the p ho tomet ric metallicity ac¬ 
curate to 0.10 dex (|Nemec et al.||2013|). 

We calculated tne metallicity of the third star, 
EPIC 210282474, to be [Fe/H] = -2.62 ± 0.15. We note 
that if the star is indeed modulated, the partial phase 
coverage of the Blazhko effect could affect the photo¬ 
metric metallicity value. We increased the uncertainty 
to account for this additional factor. Overall, the three 
metallicity values are in good agreement with the aver¬ 
age metallicity values of Leo IV (around -2.3 or -2.6) 
and confirm its very metal-poor nature. 


5. CONCLUSIONS 

We detected the variations of all three extragalactic 
RR Lyrae stars in the K2 mission Campaign 1 data. 
The stars are me mbers of the ultra-fain t dwarf spheroidal 

These are the 


2009). 


galaxy Leo IV (Moretti et al. 
faintest pulsating stars measured with the Kepler space 
telescope so far, at a brightness level of Kp « 21.5 mag. 

One of the variables, EPIC 210282473, displays clear 
amplitude variations with a period of 29.8 ± 0.9 days, 
while EPIC 210282474 is likely modulated with a period 
longer than the campaign length (80 d). In the case of 
the former star, almost three modulation cycles are cov¬ 
ered, during which the strength of the Blazhko effect and 
the relation between the amplitude and phase modula¬ 
tion both changed. These observations are the farthest 
measurements of the detailed parameters of the Blazkho 
effect at a distance of 154 kpc, and the first one beyond 
the Milky Way and the Magellanic Clouds. We deter¬ 
mined the photometric [Fe/H] indices of the stars that 
range between -2.16 and -2.64, in agreement with the 
spectroscopic values measured in Leo IV. 

We also demonstrated that with fine-tuned algorithms 
based on image subtraction, nearly the same faintness 
level can be achieved for stationary targets that is attain¬ 
able for moving objects, such as TNOs. We found that 
unlike brighter sources, very faint but high amplitude 
variations do not significantly contain frequencies caused 
by the attitude correction maneuvers of Kepler. The im¬ 
plications of this exercise are far-reaching, as photomet¬ 
ric space missions will continue to stare at crowded fields 
in the future. Globular clusters have already been mea¬ 
sured during the K2 mission, and Kepler may observe 
Cepheid stars in another galaxy, IC 1613, in Campaign 
8. The TESS and PLATO missions will employ even 
larger pixels that Kepler did, leading to strong crowd¬ 
ing in several fields wh ere image subtract i on will be a 
neces sary requirement (Ricker et al. 2014 Rauer et al. 
2014). Just one, important example is the Large Mag¬ 
ellanic Cloud within the southern continuous observing 
zone of TESS. The LMC contains several bright, large- 
amplitude variables (Cepheids for example) that will be 
accessible for that mission, but their photometry will 
present challenges very similar to the K2 observations 
presented here in terms of crowding and signal-to-noise 
ratios. 
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